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Abstract. We investigate operator splitting methods for a special class of 
nonlinear partial differential equations with delay. Using results from the the- 
ory of nonlinear contraction semigroups in Hilbert spaces, we explain the con- 
vergence of the splitting procedure. The order of the convergence is also given 
in some important, linear special cases. 



1. Introduction 

Partial differential equations with delay play an increasingly important role in 
modeling physical, chemical, economical, etc. phenomena, since it is quite natural 
to assume that past occurrences effect the model. For further motivation see for 
example the monographs by Wu [26] or Batkai and Piazzera [3]. 

There has been lots of work describing the asymptotic behavior and regularity of 
solutions, as well as in the numerical analysis of ordinary differential equations with 
delay, see for example the monograph by Bellen and Zennaro [4]. The numerical 
analysis of partial differential equations with delay, however, seems to be in its 
infancy. The present paper aims to contribute to this topic in analyzing an operator 
splitting procedure for nonlinear partial differential equations with delay. 

Operator splitting is widely used in numerical analysis of complex systems. The 
idea is to decompose the differential equation into simpler equations which can be 
solved in an effective way, and then represent the solution of the original equa- 
tion using product formulae. For ordinary differential equations, the theory seems 
to be quite complete, as witnessed in Hairer-Lubich- Wanner [15, Section 11.4,5]. 
There has also been enormous progress in the theoretical investigation of splitting 
procedures for infinite dimensional systems in recent years, see for example the 
monographs by Farago and Havasi [13], Holdcn ct al. [17], Hundsdorfer and Verwer 
[18], or Lubich [24]. See also the recent papers by Batkai et al. [2, 1], where nonau- 
tonomous equations and spatial approximations are also considered. Unfortunately, 
abstract results analyzing the order of convergence are rather incomplete, and, as 
it seems, can be applied to delay equations only with enormous difficulty. 

The idea to apply splitting procedures to delay equations is the following. Con- 
sider, e.g., a delay equation of the form 



u'(t) = Au(t) + sin \u(t - 1) + / a ■ u(t + a) da 
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with some initial and boundary conditions. The delay term appearing here repre- 
sents the two main classes of possible delays in applications: point delays corre- 
sponding to dependence on a single event in the past, and distributed delays (given 
by an integral term) corresponding to dependence on a whole time period in the 
past. In our opinion, distributed delays arc more realistic in modeling. 

Since the delay term is in a way a scalar operator, it is natural to decompose the 
equation into two sub-problems: the heat equation 

w'(t) = Atu(t), 

and a scalar-valued delay equation 

v'(t) = sin I v(t - 1) + / a ■ v(t + a) da 



-0.5 

Both equations can be solved numerically in an effective way. Note that the second 
equations becomes here an ordinary differential equation with delay, hence the 
methods described in Bellen and Zennaro [4] can be applied. 

Our aim is to put this example in an abstract perspective explaining the con- 
vergence and analyzing the order of convergence in some special cases. To be more 
specific, we will consider on the Hilbert space H the following abstract delay dif- 
ferential equation 

^- = Bu(t) + $u u t>0, 

u(0) = x e H, 

u o = feV([-l,0];H), 

where u : [—1, oo) —> H is unknown and the history function ut defined by u t {a) := 
u(t + a) for s £ [—1,0]. The further precise definitions and assumptions will be 
made in Section 2, see also the monograph Batkai and Piazzera [3]. In many cases, 
as explained above, it is easier to solve the equation without delay and the "pure" 
delay equation separately. In this case, it is natural to apply some operator splitting 
procedure described below. 

Let us fix a time-step h > and solve first the equation 



dt 

( SP1 A) \ vM(0)=x=:z 1 

v^=f=:f, 
Then for y\ := v^(h) we solve the equation 



H\ te[o,h], 



(SPl/2) 

To initialize the next step we set 



dt 

w {1 \0) =:yi. 



t& [0,h], 



x 2 := w w (h) 
f (i) 

Then we start the procedure again by solving Equation (SP1) on the time-interval 
t 6 [h, 2h] with initial values (^) . We iterate this procedure and solve the purely 
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delayed eqaution 

dwW W te[{k-i)h,kh], 



(SPk/l) 



di 

(k) , 
v {k _ 1)h = fk- 



Then we take y k := v^ k \kh) and solve equation without delay 

(SPk/2) < ^ te[(*-i)M*], 

«;W((*-l)ft)=:i/fc. 
Finally, we initialize 

x k+ i := w {k \kh) 

/fe+l : = W/cft ■ 

The sequentially split solution at time level t = kh is then 

u sq (kh) := Xfc+i. 

We shall show that for fixed i G [0, oo) and for h — > (t = kh, so k — >■ oo) this split 
solution u sq (kh) converges to u(t). 

This procedure is especially useful, if we can drastically reduce the computational 
complexity of the problem by the splitting. This is, e.g., the case if 4> = C6(— 1) is 
a point delay. Hence, we can integrate the first split equation explicitly, reducing 
the problem to solving the second equation, a classical partial differential equation. 
For a different splitting rocedure, designed specifically for distributed delays, we 
refer to Csomos and Nickel [8]. 

Besides the sequential splitting there are other splitting methods known in the 
literature. However for the sake of simplicity, in this paper we only deal with the 
sequential and the first order Lie splitting. The convergence of higher order splitting 
schemes can be explained along the same lines. More refined convergence analysis 
of higher order splitting, and numerical experiments are subject to further research. 

In Section 2, we show a way to rewrite a delay equation as an abstract Cauchy 
problem, and show how to associate nonlinear contraction semigroups to that. Also 
basic facts and about such semigroups are recalled there. These results are then 
used to explain the convergence of the splitting procedure. In the final Section 3 
we study the order of convergence in some important special, linear cases. 

2. The delay semigroup 

Consider the abstract delay equation of the following form (see, e.g., Batkai and 
Piazzera [3]): 

^ = Bu(t) + t > 0, 

W \ u(o) = x e H, 

on the Hilbert space H where I < p < oo and Ut is the history function, i.e. 
u t (a) = u(t + a), t > 0, a G [-1,0]. 
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It is possible to transform the delay equation (1) into an abstract Cauchy prob- 
lem, see Batkai and Piazzera [3]. In order to do so, we take the product space 
£ := H x L p ([— 1, 0]; if) and the new unknown function as 

t^U(t) := ( Uit) )e£. 



Then (1) can be written as an abstract Cauchy problem on the space £ 

{ U(0) = (*)e£, 

where the operator Q is given by the operator matrix 

(3) ^ 

on the domain 

D(g):={(J)efl(B)xW^([-l,0];l); /(0) = x] . 

It is shown in Batkai and Piazzera [3, Corollary 3.5, Proposition 3.9] that the delay 
equation (1) and the abstract Cauchy problem (2) are equivalent, i.e., they have 
the same solutions. More precisely, the first coordinate of the solution of (2) always 
solves (1). 

Now the abstract Cauchy problem (2) can be solved by semigroup theoretic 
methods. These we recall briefly 

The best Lipschitz constant of a Lipschitz continuous function T : X — > X on a 
Banach space X is denoted by ||T||Lip- The operator T is called then a contraction if 
1 1 r 1 1 Lip < 1. A subset A C X x X is called dissipative if for each a > 0, (I — aA) -1 
is a function and a contraction. We call A m- dissipative if D(A) = X and for 
all a > ran(i — aA) = X (D(A) is the domain of definition of A). A family 
T of Lipschitz continuous operators on the Banach space X is called a strongly 
continuous semigroup of type ui £ K if 

(1) T(t) : X X for alU > 0, 

(2) T{t)T{s)x = T(t + s)x for all t, s > 0, and T(0) = I, 

(3) ||T(i)||Lip < foralU>0, 

(4) limT(t)x = x. 

The central point in this theory is the theorem of Crandall and Liggett [7] general- 
izing the classical Hille-Yosida Theorem. 

Theorem 2.1. Let A C X x X be cu-m- dissipative for some lo £ K, i.e. suppose 
that A — cuI is m- dissipative. Then A generates a strongly continuous semigroup of 
type uj given by the formula 

(4) T(t)x= lim (I - ±AY n x 
for all x £ X and t > 0. 

If A generates a semigroup in this way, we call the abstract Cauchy problem 
associated to A well-posed on X. We return to our Cauchy problem (2), which 
under the next conditions becomes well-posed on £ . 

Assumption 2.2. Suppose that 

(1) there is a 6 R such that B— al is m-dissipative and hence B is the generator 
of a strongly continuous (nonlinear) semigroup V on H, and 
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(2) $ : W 1 ' P ([-1, 0];H) -> H is an operator defined as follows. Let g : H ->■ # 
be a Lipschitz continuous function with Lipschitz constant /3, and let 77 : 
[— f , 0] — > ^f(H) be a function of bounded variation with values as bounded 
linear operators. Assume further that 

= and that lim r,(a) ^ 0. 

er— >— 1 

Now the operator $ is given as 

$(/) := .9 f ^ d V (a)f(a) j for / e C([-l, 0];#). 

Following Webb [25, Section 4], it is possible to find a new equivalent norm in 
£ such that the operator Q — jl becomes an m-dissipativc operator and hence the 
generator of a (nonlinear) semigroup. We define 

\d V (a)\ 

-1 

denote the total variation of 77 on [— 1, r]. Let us introduce the new equivalent norm 

as 

1 



\ p 

W + J\f(<r)\ p T(a)da) for (*) E £ . 



Then we have the following result, sec Webb [25, Proposition 4.1]. 

Theorem 2.3. If p > 1, then Q — 7Z is m-dissipative in £ for 

7 = max{0,r(0)(i + f )+a}, 

where | + i = 1. If p = 1 and (3 < 1, then Q — jl is m-dissipative in £, where 
7 = max{0, t(0) + a}. 

Let us turn our attention to the splitting procedure described in the introduction. 
On the semigroup level in the product space £ this corresponds to the splitting 

. (B 0\ . _ fO $ 

with 

D{Ai) = D(B) x L p ([—1,0],t;H) 

and 

D(A 2 ) = { (*) G H X W^([-l, 0], r; ) : /(0) = *} . 
These operators generate strongly continuous semigroups of the form 

*(.)-("?>?). r,,„C)-(t» 

with / being the identity operator on L p ([— 1, 0], r; i?) and 

i 



= x + y da, 




see Webb [25, Proposition 5.12]. If we denote the projection to the first coordinate 
with 



() 
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then the sequential splitting (after k steps with time step h) can be written as 

u s *{hk)=V l [{T l {h)T 2 {h)) k { x f ) 

Our first aim is to show the convergence u sq (hk) — > u(t) = ViU(t) for k — > 00 and 
kh = t. 

The main abstract technical tool in investigating splitting procedures are Lax- 
Chernoff type theorems and variants of the Lie- Trotter product formula for non- 
linear semigroups. Such result were proved in the paper by Brezis and Pazy [5, 
Section 3], and improved by Kobayashi [20]. We recall here the important results 
for our investigation. 



Theorem 2.4. (a) Let A, B, and A + B are uj-m-dissipative sets in a Banach 
space X generating the semigroups S, T , and U , respectively. Then 

_ -i _ -, -1 n 

t/ ^ a; = „ Um o[( / -« A )" \ x 

for all x G X . 

(b ) If, in addition, X is a Hilbert space and A + B is closed, then 

u ^ = ^js{i)T{^)] n x. 

We shall refer to the case (a) as the Lie- splitting, and to case (b) as the sequential 
splitting. For more general product formulae we refer to the seminal paper by Lions 
and Mercier [23]. 

The extra condition about the Banach space in (b) is important, as a coun- 
terexample by Kurtz and Pierre [22] shows. Nevertheless, one can fairly relax 
this condition by requiring that X has uniformly Gateaux differentiable norm, see 
Kobayashi [21]. It seems to be folklore that a separable Banach space can be equiv- 
alently renormed such that the new norm has this differentiability property. Note 
however, that for a Hilbert space H and 1 < p < 00, the norm of the ii-valued 
Bochner space L P (H) is uniformly Gateaux differentiable (by a result of Day [9] 
the space L P (H) space is uniformly convex and uniformly smooth, which is even 
more than what is required; see also Diestcl [10, Sec. 2. 3]). All in all we note that 
the above theorem is applicable in the situation of the delay equation, i.e., for the 
semigroups 7i and T2 on the product space £ = H x L p ([— 1, 0]; H), as described 
above. Notice too that one has to endow the product space with an appropriate 
smooth and uniformly convex product norm (see Clarkson [6]). 

After all these preparations, we obtain the following general convergence result. 

Theorem 2.5. Suppose that Assumption 2.2 holds. Consider the delay equation 
(1) and the splitting of the operator Q = A\ + A2 described above. For every 
p G (1, 00) the sequential splitting given by the formulae (SPk/1) and (SPk/2) and 
Lie splitting converge in H x L p ([— 1, 0], r; H). 

Proof. By Theorem 2.3 the semigroups 7i, T2 and T are all of type 7 for some 
7 G K. Since Ai + A2 is readily seen to be closed (and generates T), Theorem 2.4 
and the paragraph thereafter yield the proof. 

□ 

3. Splitting in linear case 

We now turn to the case of linear delay equations and prove convergence of 
the sequential splitting together with error estimates. Further analysis of more 
complicated equations is subject to ongoing research. For the terminology on linear 
semigroups we refer to Engel and Nagel [11]. 
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First wc show how an abstract semigroup result can be applied directly to the 
problem, demonstrating the power of the semigroup approach. To obtain some error 
estimates one can apply a general result by Hansen and Ostermann [16, Theorem 
3.1]. 

Theorem 3.1. Assume that operators Ai and A 2 generate (linear) contraction 
Co-semigroups in the Hilbert space £ and that Q = A\ + A 2 generates the strongly 
continuous semigroup T and satisfies 



(6) 



D(g 2 ) c D{AiA 2 



Then the for sufficiently smooth initial values the Lie- and the sequential splittings 
have first order convergence, i.e., for every 6 D(Q 2 ), 



(I-hA 1 )- 1 (I-hA 2 )- 1 { x } )~T(h)^) 


<C -h 2 ■ 


(/) 


and 

G 2 




Ti(h)T 2 (h){*)-T(h){*) 


<C -h 2 ■ 


(/) 


G 2 ' 



Ax 



For other abstract error estimates, see for example Jahnke and Lubich [19] or 
Farago and Havasi [12]. 

Theorem 3.2. Assume that B generates a linear contraction Co-semigroup and 
that the linear operator $ satisfies the condition ran$ C D{B). Then for yA G 
D(Q 2 ) the sequential splitting is of first order. 

Proof. To prove the desired estimate, we have to check the conditions in Theorem 
3.1 for the operators 

r B 0\ . _ fO $ 

with their respective domains, see Section 2. 

By Theorem 2.3, we know that each operator Ai generates a contraction semi- 
group in £ . Hence, we only have to check the domain condition. 

Consider first the right-hand side. Then 

D{A^A 2 ) = {(p e D(A 2 ) : A 2 {% e £>(^)} 

= {(*)eD(A 2 ) : §/ e D(B)} = D(A 2 ) 

since ran$ C D(B). 

For the left-hand side wc obtain 



( B *+* f ) € D(G)} 



D(G 2 ) = D{9) : G( x f ) € D{Q)} = {(;) G D{Q) : 

{ (p g D{B) x W 2 *([-l, 0],r; H) : /(0) = x, /'(0) = Bx + */ G D{B)) . 
Hence, 



D(Q') C D(Q) c D(A 2 ) = D(AiAa), 



and the assertion is proved. 



□ 



We now want to weaken the assumption about ran3> and still obtain a conver- 
gence rate result. To do so we directly calculate the local error of the splitting. The 
splitting procedure is certainly stable, as all appearing semigroups are of type 7. 
So without changing the stability of the splitting, we can take arbitrary equivalent 
norms on £, and so get rid of the function r. Let us record this fact for later 
reference: there exist constnats M > 1 and wet such that 

(7) ||(7"i(i)r 2 (t))' £ || < Me ukt holds for all t > 0, k G N, 
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where 7i and 72 are the semigroups generated by the operators 

with D(Ai) = D(B) x LP([-1,0];F) and 

D(A 3 ) = {{*) GH xW^([-l,0];lf) : /(0) = *} , 
and are give by (5). 

Further, by the calculations in Section 2 and by the results in Batkai and Piazzcra 
[3, Section 3.1], we know that the semigroup T generated by (G,D(Q)) is given by 

(8) r(t)$ - (u{t) " 

and we have 

(9) ut(a) 




Here is now the assumption replacing the one "ran$ C D(B) V . The price we 
shall pay for weakening the assumption is, however, that we will have to improve 
the regularity of the initial condition. 

Assumption 3.3. Let Y be a Banach space, p g (1, infty), and consider the linear 
operator 

$ : W hp ([-1,0];Y) -> Y 

given by 

o 



with 77 : [0, 1] — > Jzf(Y) of bounded variation, such that for some <7o < 
77 : [ao, 0] — > J£(Y) is Lipschitz continuous. 

a) Suppose that the above holds for Y = H. 

b) Suppose that the above holds for Y = D(B). 

c) Assume further that 

Tf(— 1) = and that lira r)((r) 7^ 0. 

cr— >— 1 

Informally, this means that $ maps D(i?)-valued functions into D(B) and that 
there are no point delays accumulating around zero. Note that if the operator $ is 
"scalar" , for example if $ = cS-i with c £ C, then this condition is always satisfied, 
while the conditions of the previous Theorem 3.2 are not. 

In what follows, we abbreviate the norm notation || - ||lp([-i,o];£)(b)) to |H|lp(d(b))i 
and similarly to all function spaces defined on [—1,0] with values in H or D(B). 

Theorem 3.4. Assume that p > 1, B generates a linear contraction semigroup in 
the Hilbert space H, and that Assumption 3.3 a) and b) are satisfied. Define 

2? :={(/) ei'(S > )xW 1 '([-l,0];J)(5)), / e Lip([-1, 0]; 7J), f(0) = x}. 

Then for (%) G T> the sequential splitting is of first order. More precisely, for all 
To > there is a constant C > such that 
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where for n) G T> one defines 



\x\\ B + \\Bx\\ B + \\f\\wi-P(D(B)) + ||/||Lip(ff)- 



The proof of this theorem follows the standard route of stability analysis and 
estimating the local error using the variation of constants formula. 
These are carried out in the scries of lemmas 

where wc first show the invariance of the subspace T> under the delay semigroup, 
and then show the local error estimate for initial values from this set. 

As before, denote by V\ : £ —> H and T>2 ■ £ — > L p ([— 1, 0]; if) the coordinate 
projections. 

Lemma 3.5. Let (p G D(G) be such that f G Lip([— 1, 0]; if) with Lipschitz con- 
stant L. Then for all t G [0, 1] the function u t := T > 2T(t)(j) is Lipschitz continuous 
with constant 

max{L,M||(*)|| a }, 

with M depending only on T ■ 

Proof. By (9) we know that 

Ut : [— t, 0] H, u t {a) = u(t + a). 

This is a Lipschitz continuous function, since it is continuously differcntiable, with 
constant 

sup \\ViT(t + a)g{*)\\<M\\g(f)\\. 

o-e[-t,o] 

Now, since / is assumed to be Lipschitz continuous and since /(0) = x = u t (—t) 
(for (p £ D{Q)), and also by taking into account (9) again, it follows that u t is 
Lipschitz continuous with the asserted constant. □ 

Lemma 3.6. Let ( x f ) G £, let v t := Pa (7a (t) ( ) and u t := V 2 (T(t)(^)). Then for 
all t G [0, —do] we have 

$v t = §u t + 0(t), 

where 0(t) denotes a term of magnitude t ■ Cmax{||/||LP(.ff), ||^||} with a constant 
C independent of x and f. 

Proof. Similarly to (9) we have 

\(t + a) if t + er > 0, 



v t ((j) 
We now can write 



||$M t - $V t \\ = 



f(t + a) if t + cr < 0. 



d?7(cr) (u t (a) - v t (a)) 

< J \\u t (a)-v t (a)\\\\d V (a)\\+ J \\u(t + a) - v(t + a)\\ \\d V (a)\\ 

-1 —t 
-t 

< /||/(t + <7 )-/(t + (r)||||d»7(a)||+ f\\u(t + a) - v(t + a)\\ \\d V (a) 



n 

= 0+ / ||«(t + «T)-t;(t + a)||||d»7(a)||<2Mmax{||/|Up (H) ,||!r||}t, 
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where M depends on the bounds of the semigroups T, 1~2 and on the Lipschitz 
constant of 77. □ 

Lemma 3.7. a) The subspace T> is invariant under the semigroup T. 
b ) For G D the function 

[0, -<7 ] 9SH $1i s £ 

zs Lipschitz continuous with constant C(\\x\\b + \\f\\up(H) + ll/llw^Cir))- 
c j For G I? i/ie function 

[0, -fro] 3 s h> G D(B) 

is bounded in norm, and the function 

[0, -cr ] 3S4 F(cr)$u s G i? 

zs Lipschitz continuous with constant C(\\x\\b + || /||w 1 p(I5(b)) + ||/||Lip(H) + 
\\Bx\\ B ). 

Proof, a) Let T be the semigroup generated by the part Q of Q in £ := D(B) x 
L p ([— 1, 0]; D(B)). This semigroup exists by the general facts presented in the 
beginning of Section 2 and by Assumption 3.3.b). The domain of the generator is 

D(G) = {(*) efl(B 2 )xW^(K0];D(B)), /(0)=x>. 

Hence 

D = i>(a)n(irxLip([-i > o];iO). 

The space D(£7) is invariant under T. It is easy to see that for (^J G T> we have 
T(t)(f} = f(t)( x f ) for all t > 0. This follows from the fact that V C £>(£?), from 
the uniqueness of the classical solutions of (1), and from the comparison of the 
topologies of 8 and £. Next, since D(Q) C D(Q) from Lemma 3.5 we obtain that 

T(t)( x f )ev. 

b) Let < s < t < —(To and take (^) G V. By the definition of u t we have 



\®u s - $u t \ 





dr?(cr) (u s (a) - u t (a)) 

-1 



" 

< / Hu.CaJ-UtWIIIId^) 



-1 



-1 



\f(s + a)-f(t + a)\\-\\drj(a) 
f H/Cs + ^-^ + ^H-lldr?^) 



+ / |K« + <r)-u(t + <r)|H|d77(<r) 



The first of these three terms can be estimated as 

-/ 



|/(s + a) - f(t + <t)\\. \\d V (a)\\ < ||/HLip(H)|M|BV(J*f(ff)) ' (t ~ *)■ 



Similarly, since u is a classical solution of the delay equation (by V C D(Q) C 
D(Q)), it is continuously differentiable and hence Lipschitz continuous. Hence the 
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third term can be bounded in a similar way as the first. For the Lipschitz constant 
here note that 

and hence 

L<M(\\Bx + <i>f\\ + \\f\\ wl , P{H) ). 

Finally, in the second term we integrate on the interval [— t, — s] with t,s < — <jq. 
Hence this term can be bounded by C(t — s)(||a;|| + ||/]|lp(.H")); where C depends on 
T and on the Lipschitz constant of 77 on [<ro, 0]. 

c) By the considerations in part a), we have that for (^) G T>, the function 

is continuous, and hence especially bounded on [0, 1] by 
sup ||P a f (*)(*) || w i, (J3(B ))< sup \\T(s)( x f )\\g<M(\\Bx\\ B + \\f\\w^(D(B)))- 

sG[0,l] ' s£[0,l] 

This shows the first assertion, since <f> e _Sf (W^Q-l, 0]; D(B)), D(B)) by As- 
sumption 3.3.b). The second statement is a trivial consequence now, since for large 
A > the function [0, 1] 3 s 4 V(s)R(X, B) £ ^(H) is Lipschitz, hence so is 
s 1— > V(s)R(X,B)(X — B)$u s if we take into account also part b). The Lipschitz 
constant is 



L — C[\\x\\b + \\f\\w 1 -P(D(B)) + 
where C does not depend on x and /. 



iLip(if) 



+ \\Bx + $f\ 



□ 



Lemma 3.8 (Local error). There is a constant K > such that for allt € [0, — 00) 
we have 



Ti(t)T 2 (t)(*)-T(t)(* f ) 



< Kt 2 



for all (*) e V. 



Proof. Take (^) 6 2?, and let w, -u t , v, v t be as in (5) and (8), respectively. We can 
write 

t t 
V(t)v(t) - u{t) = V{t)x + V{t) J $v s ds - V(t)x -J V(t- s)$u s ds. 



By Lemma 3.6, this is further equal to 

t t t 

= V(t)x + V(t) J 1>u s ds - V(t)x -J V(t- s)$u s ds + J 0{s)ds 



V(t - s) {V(s)$u s - $w s ) ds + 0(t 2 ) 



V(t - s) (V(s)<i>u s - $w ) ds + V(t- s) ($m - $u s ) ds + 0(t 2 ). 



By Lemma 3.7, the function [0, — <jq\ 3 s t— > <£>it s is Lipschitz-continuous, and the 
function [0, — (To] 9 s 1— > V(s)$u s is Lipschitz continuous. Denoting their Lipschitz 
constant with L, we obtain that 



\\V(t)u(t) - v(t)\\ < MLt 2 + MLt 2 + \\0(t 2 )\\ < Ct 2 
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As for the second coordinates, we can write 
o 



u t - v t \ 



LP (H) 



\\u t (a) - v t (a)\\P da 



-1 
o 



\\u(t + a)- v(t + a)\\ p da < Ct 2 ^ +1 (\\f\\ hP{H) + \\x\\y. 



This finishes the proof by defining K in terms of the constant C from the above. □ 

Proof of Theorem 3.4- Without loss of generality we may assume Tq = 1 (the con- 
stant C shall depend on To, however). By a standard telescopic summation we 
get 

(10) (Tx(i)T 2 e)" - T(t) = £ 1 (T 1 (^)T 2 (^))"- fc (Ti(£)T 2 (£) - T(±)) T(\ 



■ kt ' 



k=0 



As we have seen in Lemma 3.7, the set T> is invariant under the delay semigroup 
T, and also from this Lemma 3.7 and Lemma 3.5 it follows that for some M > 1 



T(s){ 



< M 



T< 



for all s G [0,1], (*) G V. 



From Theorem 2.5 we obtain a possibly larger M , which makes it possible to esti- 
mate the first factors on the right hand side of (10) by exponential terms. Hence, 
by applying Lemma 3.8 we obtain that 



((Ti(£)T a (£))" - T(t)) (p ^ < E T(M)(^ 
where C* = M 2 Ke^. 



k=0 



C 
< — 

T> n 



v 



□ 
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